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ABSTRACT 



Two prediction models for graphical prognosis are developed and 
evaluated. The primary feature of these models which represents a de- 
partiure from previous graphical methods is the inclusion of the ter- 
rain-induced vertical velocities. These vertic&l motions may in turn 
produce vorticity changes, and thus height changes. 

The first model presented is con5>rised of two separate partsj the 
first part is a prediction of the ^00-mb surface and the second part a 
prediction of the 1000-mb to 500-mb thickness. The two parts may be com- 
bined to produce a svirface (1000-mb) prognosis. This model failed to pro- 
duce satisfactory results, and hence a second model was developed. 

The second model presented follows somewhat the development of the 
first, but avoids an assur^jtion relating the surface and 500-rab winds 
that is believed to have contributed largely to the original failure. 

This model operates on the current 1000-mb and 1000-mb to $00-mb thick- 
ness charts; and, utilizing a fraction of the 500-mb height change as 
predicted by arc'- other established method, produces a surface (1000-mb) 
prognosis. The model was tested in two ways; one xid.th the terrain factor 
included in the prediction equation and one omitting the terrain factor. 
Both methods produced satisfactory prognoses. However, inclusion of the 
terrain factor consistently produced the better results. 

The writer wishes to express his appreciation for the assistance 
and guidance of Professor George J. Haltiner of the U. S. Naval Postgrad- 
uate School in the preparation of this paper. 
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PREFACE 



Graphical methods of rmraerical prognosis have recently received 
vri-despread attention since the development by Fjortoft £l] of a predic- 
tion model for the ^OO-mb surface. Fjortoft provided a graphical tech- 
nique for integration of the barotropic vorticity equation. Despite the 
relative simplicity of the method and the rather restrictive assumptions 
inherent in its development, tests by the Air Weather Service of the 
United States Air Force indicate accuracy conparable to that achieved 
by a trained forecaster using conventional techniques. 

Following the success of Fjortoft 's model, several other more in- 
volved graphical models have recently been devised. Although none as yet 
have been subjected to the degree of testing that Fjortoft ’s model has, 
preliminary results are encouraging. 

The graphical methods devised thus far appear to have been least 
successful in mountainous terrain areas. The reason for this probably 
lies in the simplified boundary condition, 00 = 0 at p^, which has been 
assumed in the previous models. 

It is the purpose of this paper to attempt to incorporate terrain 
effects into graphical techniques of prognosis. The development here is 
somewhat similar to the treatments of Estoque [3^ and Reed fijj . 
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CHAPTER X 



FIRST PREDICTION MODEL 



The first model is conprised of two peirts; the prognosis of the 
500-mb surface and a separate prognosis of the 1000-mb to 500-nib thick- 
ness. The two prognoses may be combined for a resultant surface (lOOO-rab) 



The first teiro on the idght represents the contribution due to the oro- 
graphic effect and is assumed to vary linearly with press^^re (p/p^), 
from a maximum at p^ to zero at the top of the atmosphere, k is a scalar 
multiplier T-iith possible values ranging between zero and unity. It is 
assumed to equal unity for this development and hence will not appear 
in subsequent equations in this chapter. The second term on the right 
is the contribution due to the non-orographic effects, and is assumed 
to vary parabolically from zero at p=rp^ to a maximum of 0.)^^ at p = p^. 
Inherent in this distribution is the assumption that p^ (to be taken as 
^OOmbs) is a level of non-divergence where there are no orographic effects. 
The non-orographic contribution assumed here varies slightly from Estoque's 
assumption of a sinusoidal variation of vertical velocity . 

Differentiating (l) vrith respect to pressure yields 



prognosis. The development 




1 . Development of the prediction equation for the 500 -nib surface. 



The "vertical velocity" is assumed to vary thusly: 







( 1 ) 




(2) 



1 



Evaluating (2) at gives 

( ^ — 

U?4" 



ft 



vi-y'H 



(3) 



The following fom of the vorticity equation is assumed; 



If +V-Vl=-f|f 

dp 



(h) 



This form of the vorticity equation neglects the vertical advection of 
vorticity, the vortex tube term, and ^ Since = ~V • V , Eq, 

' ^ p ' C7 r 

(3) shows that the only divergence considered at ^OOmbs is due to ter- 
rain-induced vertical velocities. 

Combining (3) and (U) yields 

Vs-VH . (5, 

In order to obtain conservation equations utilizing only 500 -mb 
data, the surface wind must be eliminated from Eq. (5), This will be 
done by assuming that Vg— aV, where a is a scalar multiplier. Then 
(5) may be written 



— +V-VX = - V- V H . 



(6) 



Taking mean values of f and for the area to be prognosticated, 
(6) may be written 

+y-v(^a + . 






( 7 ) 



With the assumption of geostrophic vorticity, this becomes 

|| + v-v(f VZ +-f +^ ^'^ -h)=0 ; 

at -V- — \t 1^^ / 



or 



2 




as f and agf/gH/p^ are Invariant with time. 

Introducing the finite difference foiroula for vorticity into Eq, 
( 8) gives 




M-f 






o 



-td* 

Multiplying through by factor fdV^^g> which is invariant vrith time, 
and ejqjanding, produces 



(9) 
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(a) 



(2-2} + -^- 4 - + — — H 

M Mg 4 m^0 11 

(c) 
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(b) 



(d) 



V.^ f -}- V‘ V h = o 

Considering terms of Eq, (10) separately: 

(a) This term may be written d/ 0t(Z - Z), as other terras are 
invariant with time, 

(b) The factor Iim^g/fd^ may be taken outside of V operator, as 
variation is small conpared vrith the variation of (Z - Z). 

(c) Following Fjortoft , this terra may be witten as 



V*Vd(<?)> where 



-'P 



J((f) = 



d^smf 6059 



J, 



m' 






(d) Taking mean value of f for area to be prognosticated, this 
term may be vrritten 
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v-v 



a 



T/^cJ 



H 



3 



Rewriting (10) with above changes gives 



dr 






-— 2 . — 2 . 

d 



Mm fo 



H 



- o 



With the notation 



M m to 



Eq» (11) becomes 

dX 



:;^(2-2) +V'-'7rZ-Z-l-J(<f) ^B1- 

BX ^ J 



- o , 



As B and J(<p) are invariant with time, (12) may be written 

^ (z -2 + j(<^) + 6 ) +y • V(z- 2 -I-Jof) + B) = O 



(11) 



( 12 ) 



(13) 



or, in Jacobian form, 

-l(z-Z+J(<f) + 5) + -f jfz, 2-Z+J(‘f) + B)=0^ uu) 

dr / \ / 

which is the desired prediction equation. However, in using the above 
equation, the advection would be accon^jlished with the current 500 -mb 
wind, which obviously jnay vary considerably with time, especially if a 
prognostic interval of 12 or 2li hours is enployed. In order to arrive 
at an advective field which is more conservative with time, the Jacob- 
ian identity J(X,Y) =J(X+ Y, Y) may be employed. Applying this to (lit) 
gives 

^ (z-z + Jc S') +- B) +|- J(Z + J(cf) + d, ?-Z +-J<v + («) 

Thus the advection is accomplished with a much more conservative field, 
(^+J(<f) -h B), being the only quantity that varies with time, and this 
much more slowly than Z, 

As a preliminary test of results of the model developed herein, 

u 



Fig, 1 was prepared. An initial west-east zonal flow of 50 knots was 
assumed at 500iribs. The surface wind was assumed to be equal in direc- 
tion and ^ the magnitude of the 500-mb wind (i.e.j 3=2). Eq, (15) 
was used and a prognostic chart for a time of six hours was prepared, 
using height changes as approscimated by the first term of Fjortoft's 
series ra. The resultant prognostic chart shows a trough in the lee 
of the Rocky Mountains, as would be expected under the circumstances. 

2. Testing of the prediction equation for the 500-mb surface. 

The first step in the employment of Eq. (15) was the construction 
of the 1 . chart. This chart was prepared by conventional graphical methods, 
using a grid distance d of lO^km as suggested by Estoque , Next, the 
conto\irs of J((f) corresponding to a grid distance of 10-^km were added 

graphically to the ^ chart. Assuming a value of ^ for a, the quantity B 

—2 

was evaluated as 1.5 x 10 , assviming standard values of the meteorolog- 

ical parameters, neglecting variations of m^, and taking d as lO^km, as 
above. (Variations of m are small, as a conforroal projection with stand- 
ard latitudes of 30N and 6ON was used, and prognostication was confined 
to middle latitudes) . A terrain contour chart of North America was multi- 
plied by the above constant and isolines of the quantity B were drawn. 

The B chart was then added to the Z + J(*f) chart to produce the advective 
field (^ + J('f)+ B) » The Z chart was then subtracted from the advective 
field to produce the (Z - Z + J('f) B) field, or the conserved quantity 
to be advected. 

The above two charts v;ere superimposed in register, and isolines of 
(Z - Z + J(<f) + B) advected geostrophicolly with the (Z + J(^f) -h B) chart 
for the chosen prognostic interval of 12 hours. 
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The prognostic 12 -hr height changes may have been recovered by use 
of Fjortoft's series []]3 > 

A H = A (f -Z + J«f) + B) + 2(f +Jof) + 

and compared with observed values of height changes, but for this test a 
more direct method of evaluation was used. The quantity (^ - Z + J(tf) + B) 
was advected for the 12-hr period and con^^ared vdth the observed values 
of the same quantity as computed from the observed chart at the end of the 
prognostic period. 

The results of four tests using the above techniques and assumptions 
were generally poor. Indications were that in some isolated regions the 
inclusion of the teirain factor gave better results than the use of the 
prediction equation 

+Jccf))+ ^ J )^ ? -z -hJctf))- O. (16) 

(Except for the different vertical velocity distribution assximed, this 
equation is essentially the same as Estoque's model [3]). However, these 
regions of improved results were not constant in space nor were they pre- 
dictable, and the areas of poorer results generally equalled or exceeded 
the ijiproved regions, resulting in generally poorer prognoses by inclusion 
of the terrain factor than by ignoring it. 

It was reasoned that the choice of a=r i,e,, that the surface 
wind is equal in direction and -g the magnitude of the 500 ~ni 6 wind (in 
connection with the terrain factor), may have been largely responsible 
for the failure of the method. Therefore, another technique was develop- 
ed using a equal to either -g-, zero, or -g, as determined in the follow- 
ing manner. The current and 850-mb charts were superimposed, and 

areas delineated where the ^OO-rab contours were approximately the same 
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direction as, at right angles to, or opposite in direction to the 850-nib 
contours. The B^O-inb chart was selected as the chart most nearly represent- 
ing the s\irface in the majority of the mountainous areas. Then, in pro- 
ducing the (Z + J(<f) -f B) chart, the contours of B as drawn for the first 
technique were added, considered" zero, or subtracted in the areas indicat- 
ed, The subsequent steps in this technique were as before. 

Results of using the new technique were slightly better than for the 
assumption that a is everyvrtiere ■^, However, improvements realized were in- 
sufficient to justify use of the method. A prime difficulty in applica- 

> 

tion was the determination of appropriate values of a, "which appeared to 
vary rapidly over small areas, especially in situations where closed lows 
were present at the 500-mb or 8^0-mb levels. An extensive search was made 
to locate a synoptic situation where the differences between 500-nib and 
850-rab wind fields were clear-cut, in an attempt to evaluate the technique 
objectively. Such a situation was located, and results of the testing in- 
dicated that whereas the sign of a was correctly chosen by the above 
technique, the magnitude of a appeared to vary considerably from the as- 
svimed values of zero, and As a consequence, the prognostication 
was not aopreciably better than that using Eq. (16), and indications were 
that even if the sign of a v/ere predictable, the magnitude would vary to 
such a degree that prognostications ignoring the terrain factor would be 
more reliable. 

3. Conclusions as to prognosis of the ^OO-mb surface. 

No improvement over Estoque's prediction model ^3]] is realized by 
inclusion of terrain effects on vertical motion by methods developed in 
this paper. It is believed that the major reason for the failure lies in 
the difficulty of expressing the surface "vjind in terms of the ^OO-mb "wind 

8 



in the 500 -nib conservation equation. 



Development of the prediction equation for the 1000- to 500-inb thick- 
ness. 



A detailed development of the prediction equation will not be pre- 
sented here. The development follows closely the methods of Section 1 of 
this chapter, so a listing of the major assxinqjtions and the final equa- 
tion will be shown. 

Assumptions: 

(a) The vertical velocity distribution is as described in Eq. (1), 

(b) Temperature changes are adiabatic. The following form of the 
adiabatic thermal equation is employed: 

+COO- =0 („) 



(c) Vr»VdL ^"‘dVo'V-f are relatively small and are neglected. 

(d) Similar assviraptions as those used in progressing from Eq. (10) 
to Eq. (11) in Section 1 of this chapter, 

V/ith the notation 






and 




+ 



0 ” 



the prediction equation evolves as 

^ ('h - Ki^h + G) -f-V* y (^h 'Kqh -I- Jccf) + (i8) 



or^in Jacobian form^ 

(h -K,h+ J«f) + G) + 1 J(?/W - K,h + +C-) = 0. 

Due to the failxire of the 500-mb model to produce an acceptable prog- 
nosis, no attempt was made to test Eq. (19). Rather, a different develop- 
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ment was attempted in 'vdiich the assumption a^ was not necessary, as 

it is believed that this assumption was the major cause of failure. The 
new development is presented in Chap. XT. 



/ 



10 



CHAPTER nr 



SECOND PREDICTION MODEL 

This model is designed to be applied in conjunction with a separate 
500 -mb prognosis, such prognosis to be made by any method the user's ex- 
perience deems the most accurate. Fjortoft's method [l] is suggested* The 
development herein parallels that of Reed 



1. Development of the prediction equation. 

The vertical velocity distribution of Eq. (1) is also assumed for 
this model. Differentiating (1) with respect to- pressure yields 






d p % 

k 

and evaluation at lOOOmbs gives 



( 20 ) 



dpi % a 



( 21 ) 



Assxxming that temperature changes are adiabatic, the following form 
of the adiabatic thermal equation is employed: 



ar.[dpj^-l S 



hi. 

a p 






O 



(22) 



By integration from p to p^, (22) becomes 

U 

Ca) dp , (23) 

The surface wind may be used for the advection in lieu of a suitable 
mean wind for the layer by considering the wind hodograph to be linear 
and the thermal wind relationship to hold, i.e. : 



ah 



at 






+ Vs‘Vh=- cr 
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/ 



/ 



/ 



/ 



/ 



/ 




Evaluating the right side of (23) by use of (1) leads to 

cr^cod,. =-^Kg^(\/5-yh'):^ j °dp 



-f crco 



4-(UIh ■ 



n 



Performing indicated integration, and noting that P^— 2p^, gives 

cod^ = -o-Kg/SCys-VH) 



(T 



T 



+ 0-OJ„ 



^ 3(f,-p5) 



IV 

Zfg 



-p 

rj cudp = -^ 0- Kg (\4 • 7 H)fj + A <r P5 aJ„ . (2U) 

Jn. ' 



which finally reduces to 

% 

cr' 

% 

Substituting from Eq. (21;) into Eq, (23) produces 

'Y +|- O' '"'3 (7 - V 4 cTf, 



3't »Wv L| 

Evaluating Eq, (1;) at Pg yields 

a: 
a 



5 



^4_V5 <7rj’ +-f) - i 

at ^ Wot-Ty - t ) 



(25) 



(26) 
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assuming that = V^, This assumption is used throughout this develop- 
ment. Combining (26) and (21) produces 



10 ' ^ 



(27) 



r\ 

Multiplying (2^) by factor 3f/crp^ gives 

2'fco„ i _ 3~f 

0- vi 



g+y-Yh+io-K3/2(y-yHji| 

which, upon substituing in (27)^ yields 



( 28 ) 



■n 






_ 3-f 



dh 



3 1 



+y^'Yh+|crK9^(ys-yHjl| 



(29) 



Collecting advective terms and time derivatives, and again noting that 
Pq— 2p^, (29) may be written 

Making the geostrophic approximation for leads to 






(30) 



A 

dt 



H 






o-f; 



-y-v 



Mhp ^ 



(?o-Zo)+J^zh+-f + 

*5 






IP. 



H 



(31) 



which, after multiplying by factor fdVlh>i^g> yields 
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at 



2 ,2 



(2o h 



V,-7 



(Zo-?.) + 



3-f'd" ,^T .IIK+dVo 

h +J(cf) + 






?o 



H 



( 32 ) 



where 






j( 



V) 



o 



_n d sin cos ^ 

m^g 



df 



following Fjortoft ([iQ , With the notation 

S 

(32) may be written 



z \Z 



F= 



and 



r 3~^ d 

L ~ , z _ ^ 402. 



Hm g 






at 






tF| = Vs ' V ^o“?o+CF-tJ(c|)) + F 



as J(^p) and F are invariant with tinE. This is the prediction equation, 
which may be written in Jacobian form as 



H„-?o+Ch-(Jc<f) +l^ = - _,?yZ„tCh+J(f)+F 



(33) 



However, if (33) were to be used in prognosis, the ad’/ection would 
be acconplished vxith the current 1000-rab chart, •vrtiich may vary consider- 
ably over a prognostic period of 12 or 2li hours. In order to arrive at a 
more conservative advective field, the Jacobian identity 

J(x.Y) = J(x+Y,r) 

may be emplqyed. Applying this to (33) leads to 

a 



at 



+Jccf)+pj = 



- -f- J + F, *^Vc^)+ F 

T h 



(3U) 



lii 



which accomplishes the advection vdth a more stable field than Z^, 

To recover height changes, the local change in the advected quantity 
over the prognostic period may be written 



A 



+ Ch -f- J(cf) + F 



-■-f\ 



or 



A 



Zo'?» 







as F and J(if) are invariant with time. With the relationship h = Z 
(35) niay be written 



(35) 

5 “ 2o> 



-C 2 



= -A 



Solving for A Zq yields 



A2^ = — A -t- A2c +C AZs- 






Utilizing Fjortoft's method [^l] to evaluate gives 









C-H 



A -V CA^. 



Substituting ( 38 ) into (37) yields 



r 



A2„ = 



C-f* ' 



A A 

C-H 



+ 

C-H ^ 



Testing indicated that A-;= 2J, so (39) is simplified to 

r . 



A2. = A 






Cc-^0 



c^\ ^ 



which fiirther reduces to 



(36) 



(37) 



(38) 



(39) 



(UO) 



4 - o. b (ia) 
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For simplicity of application, Eq, (Ul) was further simplified to 

A = A + 0.SAZ5- ^ (^2) 

which was the equation utilized. However, testing indicated that Eq. 

( )[1 ) may have given better results than the simplified form of Eq. (lj.2). 



2. Testing of the prediction equation. 

The prognostic period of 12 hours was chosen to test the prediction 
equation. Four successive 12-hr intervals v;ere prognosticated starting 
with the OO 3 OZ chart of lU November, 19$$ • Since the purpose of this 
paper is to attempt inrorovement in graphical prognosis by the inclusion 
of terrain effects, the prediction equation was aoplied in two different 
forms. First, a prognosis was made applying Eq. (3U)^ and a second prog- 
nosis using 












‘f) 



-f- ^ 



r 






? 0 c h + J ( ^ 2 o + ^ ^ + J (cf: 



U3) 



which is singly (3U) neglecting vertical motion induced by terrain. With 
minor differences, Eq. (U3) is the prediction equation derived by Reed 
In appljdng (3h)t the folloi-rlng steps were carried out; 

(a) The current Yq chart was constructed by conventional graphical 
methods, using a grid distance of 600Kra. 

(b) The current thickness chart was multiplied by the constant C 
and added to the "2^ chart. In practice, using a value of 0.75 for C in 
this step was found to give better results than the derived value of 0.9, 
and consequently this value was used. 

(c) Considering k = 1, the coefficient of H in the quantity F was 
evaluated to be 7.35 x 10 A terrain contour chart was multiplied by 
this factor to produce the quantity F, and isolines of F were added to 
contours of J((f). The chart of F -+ J(if) (Fig. 2) was added to contours 
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produced in step (b) to produce the advective field ( Zq - h Ch J(tp) + F), 

(d) The current Z^ chart was graphically subtracted from the ad- 
vective field to produce the conserved quantity (Zq - Z^ + Ch + J(^) -h F), 

(e) Charts (c) and (d) were superimposed, and isolines of the con- 
served quantity advected geostrophically with the advective field for the 
prognostic period of 12 hours. A difficulty arose in the advection process 
along the west coast of the United States and Canada due to the extreme 
gradients of the advective field caused by the high coastal mountain 
ranges. It was reasoned that where the terrain gradient was very large 
the horizontal wind would not be fully translated into upslope motion, 
and therefore result in lesser vertical velocities. Thus the scalar multi- 
plier k in Eq. (32) may be appreciably less than unity in regions of such 
strong terrain gradients. The amoxint of reduction of k was determined by 
cases tested. This reduction factor may be taken into account by appro- 
priately reducing the gradients of F where the gradient is very large. 

The final form of F is given in Fig. 2, 

(f) The advected isolines of (Z^ - Zq - b Ch+ J(<f) + F) were sub- 
tracted graphically from the initial isolines of the same quantity to 
produce isolines of A, 

(g) Half of the 12-hr 500-rab height change was added graphically 
to A to produce the ^Zq chart. In practice, any suitable method of prog- 
nosticating the 500-mb chart could be employed, but for this test the 
actual 300 -mb height change was used in order that the prediction equa- 
tions, Eqs, (3U) and (Ii3)> might be more objectively compared. 

In applying Eq. (U3), the steps were identical to the above, with 
the exception that in step (c) only the contours of J(if) were added. 

The results of the above tests are displayed in Table 1, which cora- 
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TABLE 1 





1 


2 


2 




Correlation coefficient ; A 


0.832 


0.8U7 


0.710 


0.790 


B 


0.707 


0.673 


0.580 


0.630 


C 


O.SlUi 


0.851 


0.790 


0.812 


D 


0,866 


0.870 


0.92ii 


0.8U8 


Root-raean-square error : A 


9Uft 


116ft 


103ft 


99ft 


B 


ll5ft 


172ft 


169ft 


l60ft 


C 


85ft 


102ft 


88ft 


89 ft 


D 


67ft 


7lift 


53ft 


69ft 



A : Prognostication including terrain factorj couputations for 
area west of 96W, i.e., Rocky Mountain area, 

B : Prognostication excluding terrain factorj computations for 
area west of 96W, 

G : Prognostication including terrain factorj computations for 
entire area covered by gid.d, 

D : Prognostication including terrain factorj computations for 
area east of 96W, 

1 : 12-hr prognostic period commencing 0030Z lU November, 1955 * 

2 ; " " " " 1230Z " " " . 

3 : " " II n 0030Z 15 " ” . 

I " II II II I23OZ " " " , 
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pares height changes as prognosticated by the two methods to the observed 
height changes. The statistical parameters were congjuted using a ij ,8 point 
square grid, which enclosed the area bounded by points 28N, 12l|.V/j 55N, 
136W; 53N, 65Wj and 271'J, 79W, Figs, 3, and $ illustrate height changes 
for the period 1230Z 15 November to 0030Z 16 November, 1955 • Fig, 3 shows 
observed height changes. Fig, I 4 . the height changes as prognosticated by 
Eq, (3U)j and Fig. 5 the height changes as prognosticated by Eq. (h3). 

The grid encompasses an area slightly larger than the above illustrations, 

3 . Conclusions, 

The inclusion of the terrain factor in this model results in materi- 
ally improved prognoses. 

However, the process of making a prognosis by this method is rather 
lengthy, and whereas the inclusion of the terrain factor into the pre- 
diction equation adds no more steps to the routine, (F may be combined 
with J((f) on a permanent chart and used for each prognosis), the advec- 
tive procedure is nnich more delicate. The advective field has areas of 
very strong gradient mth the isolines of the conserved quantity closely 
packed and almost parallel to the advective field. Consequently consider- 
able care must be taken to’ insure accurate advection. 

In constructing the contours of F, terrain irregularities were 
smoothed to some extent. Some smoothing is mandatory, as a minutely 
accurate F chart is not reasonable to use when considering large scale 
flov7. Also, it is evident that the smoother the F chart, the easier the 
graphical constructions will be. Further testing might indicate an op- 
timum balance between detail of the terrain features and representation 
for large scale flow. 

As noted in Sect. 2, a value of 0.75 was used as a multiplier of h 
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Fig. k 
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Fig. 5 
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in producing the advective field rather than the derive^ value of 0.9. 
More extensive testing would doubtless indicate the optimum value of 
this multiplier. Reed fU] found that a value of 0.5 gave good results 
in the three cases tested with his model. 
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